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ABSTRACT 

The amount of mass in small, dark matter clumps within galaxies (substructure) is an important test of cold dark matter. One approach 
to measuring the substructure mass fraction is to analyze the fluxes of images that have been strongly lensed by a galaxy. Flux ratios 
between images that are anomalous with respect to smooth (no substructure) models have previously suggested that there is a greater 
amount of substructure than found in dark matter simulations. One measure of anomalous flux ratios is parity dependence - that the 
fluxes of different images of a source are perturbed differently. In this paper, we discuss parity dependence as a probe of dark matter 
substructure. We find that reproducing the observed parity dependence requires a significant alignment between concentrated dark 
matter clumps and images. The results may imply a larger fraction of mass in substructures than suggested by some dark matter 
simulations and that the observed parity dependence is unlikely to be reproduced by luminous satellites of lens galaxies. 

Key words, galaxies: halos - gravitational lensing - theory: dark matter 



^ ■ 1. Introduction 



, Strong gravitational lenses are important probes of the mass dis- 

■ tribution in the universe and have been used to measure both cos- 
\Q ' mological parameters and to understand and constrain the struc- 
C*~) . ture of dark matter halos and subhalos. One particular problem 

1 in cold dark matter (CDM) that may be addressed with strong 
gravitational lensing is the discrepancy between the relatively 
. small number of observed satellite galaxies around galaxies such 
t-H ' as our own and the numerous small-mass clumps seen in sim- 
OO I ulations of dark matte r halos ([R lvpin et al. 1999; Moore et al.l 
O ■ 119991: iGao et all 120041) . iDalal & Kochanek! (120021) attempted to 
^1 ' constrain the fraction of mass in dark matter substructures using 
• i— i , radio observations of quasar sources that have been quadruply- 

■ lensed by galaxies. In such systems, dark matter clumps may 
?H \ perturb the magnification of images, inducing both flux ratios 

. . . . among the images that are not well-fit by smooth lens mod- 
els and small astrometric perturbations to the im age positions 
dMao & Schneiderll998HMetcalf & Madaull200ll) . Radio obser- 
vations are necessary in order to avoid conta mination from per- 
turbat ions from stars in the lensing galaxies. IDalal & Koch anek 
(120021) found a surface mass fraction between 0.6% and 7% 
within 90 percent confidence intervals. This result is marginally 
consistent with the average surface mass density in substructures 
at the projected separation of the images as estimated by dark 
matter simulations at z = of / sub = 0.5% dMao et al.ll2004l) . 

An alternate method to investigating the substructure mass 
fraction compares observations and simulations directly, bypass- 
ing lens modeling. For lens systems with four images of the 
source, a lens configuration where three of the images are found 
near one another is referred to as a cusp. For a smooth lens halo, 
the cusp relation suggests that the sum of the fluxes in the outer 
two nearby images should be similar to that of the central im- 
age. Several papers compare deviations from the cusp relation in 
observations and simulations and suggest that there is more sub- 



structure in the obs ervations than can be accounted for in dark 
matter simulations dMaccio et al . 20061 lAmara et~alll2006l) . 

In anomalous flux ratio lens systems, the distribution of mag- 
nification perturbations is not the same for each image of a four- 
image lens: the amount by which the image is brightened or 
dimmed relative to the best-fit lens model is dependent on the 
parity of the image. Image parity describes the relative orienta- 
tion of an image to the source, and - as it is equivalent to the 
signs of the eigenvalues of the magnification matrix - all images 
can be described as having either positive parity (double-positive 
or double-negative eigenvalues) or negative parity. Parity depen- 
dence in radio observations has been used previously as evidence 
that anomalous flux ratios are due to dark matter substructure 
and not to any propagati on effect, whose effects w ould be parity 
independent in general. Kochane k & Dalall d2004l) find that the 
parity dependence in the flux anomalies is such that the bright- 
est negative parity image is preferentially made dimmer by sub- 
structures when compared to the best-fit smooth macromodel. 

On the theoretical side. iKeetonl d2001b) makes simple esti- 
mates of the effects of perturbers, ignoring the effects of pertur- 
bations to the image positions and using a linear analysis for the 
magnification perturbations. He finds that for singular isother- 
mal sphere (SIS) perturbers, double-positive parity images may 
only be brightened while double-negative parity images may 
only be dimmed and negative parity images may both be bright- 
ened or dimmed. It is possible that th is may result in an average d 
behavior similar to that observed by Kochanek & Dalai (2004). 
The parity dependence descri bed by IKeetonl d20 01b). however, 
is limited to untruncated SISs. lRozo et al.l(l2006l) . using a similar 
analysis, calculate the magnification perturbations for more real- 
istic clumps and study their behavior in statistical samples. They 
find that for a substructure mass fraction of 0.5%, the average 
magnification perturbation is small (< 1 %), and it is very slightly 
more likely that negative parity images are made brighter by sub- 
structure and positive parity images are made dimmer, an aver- 
aged result in contrast to the observations. 
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The difficulties in reproducing the observed parity depen- 
dence suggest that these simple linear perturbation estimates are 
insufficient, and astrometric perturbations are significant and/or 
substructure causes non-linear perturbations, both of which 
would require lens modeling to investigate]. In this paper, we 
attempt to reproduce the observed parity dependence by model- 
ing mock lens systems with various distributions of substructure. 
In addition to resolving the contradictions between observations 
and theoretical estimates, we attempt to constrain the kinds of 
substructure models that will reproduce the observed parity de- 
pendence. We show that the distribution of substruc t ure ne ces- 
sary is inconsistent both with the results of lMao et al.l d2004l) and 
with the anomalous flux ratios being solely due to observed lu- 
minous satellite galaxies. 

The paper is organized as follows. Section [2] describes the 
lens modeling, while Section[3]discusses the observational sam- 
ple and reproduces the previously observed parity dependence. 
Section|4]describes the substructure models used and how mock 
observations are created. Section [5] presents the results of test- 
ing several substructure models and compares them to the ob- 
servational sample. Discussion and conclusions are presented in 
Sections [6] and[7] 



2. Lens Modeling 

We model the positions of lensed images using an automated 
lens fitting code with a singular isothermal ellipsoid (SIE) mass 
distribution and an external shear component. The 9 model pa- 
rameters are 1.) the Einstein radius b; 2.) the projected axis ratio 
q; 3.) the orientation of the halo ellipticity q ; 4.) the external 
shear y which describes the effect of structure near the lens halo 
(e.g., a nearby group of galaxies); 5.) the orientation of the shear 
Q y ; 6.) & 7.) the source position, x source and Source; and 8.) & 
9.) the center of lens, x center and Renter- The SIE has a projected 
surface density 



2f 



(1) 



+ y 2 /q 2 



where £ is the elliptical coordinate satisfying £ : 
and where x and y are the Cartesian coordinates and q is the axis 
ratio. In the limit of circular symmetry, b, the Einstein radius, is 
related to the 1-d velocity dispersion <x by 



b = 4n 



lay As_ 
lei Da, 



(2) 



where c is the speed of light in a vacuum, and Di s and D os are the 
angular diameter distances from the lens to the source and from 
the observer to the source, respectively. 

We restrict ourselves to lensed systems with four observed 
images. The observed constraints consist of 8 coordinates of the 
images and a constraint from the observed center of the lens po- 
tential (presumably the position of the lensing galaxy), x cen t e r 
and Renter, for a total of 10 constraints. Fluxes are not included 
as constraints. 

The lens modeling algorithm requires several steps and is 
based upon the publicly available gravlens lens modeling code 
dKeetonll2001al) . We employ a downhill simplex minimization 
routine, and we first find best-fit parameters first in the source 



1 While the estimates of the size of astrometric perturbations are rela- 
tively small (milliarcsecond scale) dChen et alj20 07: Metcalf & Madau 
120011) . they may affect the best-fit macromodel and the magnification 
perturbations in significant ways. 



plane, which is faster than fitting in the image plane0 We start 
by finding appropriate values for the lens halo parameters - b, q, 
and q - then fixing the lens halo parameters and finding appro- 
priate values for the external shear y and 9 y . We do several iter- 
ations of fits for all the parameters in the source plane. At each 
iteration, we slightly perturb the image positions. Perturbing the 
image positions gives us several sets of data with formally iden- 
tical observational constraints and may help to avoid ending in 
local minima. Finally, we do several iterations of fitting all the 
parameters in the image plane. Additional details are found in 
AppendixlAl 

3. Observed Lens Systems 

We use a set of observed lens systems both to test the ob- 
served par ity dependence an d as t he basis for mock ob- 
servations. IDalal & Kochanekl d2002l) use 7 four-image lens 
systems with anomalous flux ratios; their sample also rep- 
resents a fair fraction of all 4-image systems with g alaxy 
lenses observed in the radio. iKochanek & Dalai (2004) use 
those 7 systems with the exception of PG1 115+080 and add 
two more systems in order to test the parity dependence of 
anomalous flux ratio sy stems. In this pape r, we u se 5 of the 
7 lens systems used in iDalal & K ochanekl d2002l) . excluding 
PG1 115+080 which does not have observed radio fluxes and 
B 1608 +656 whi ch is a double lens system. The sy stems used are 
MG0414+0534 dKatz et al] [l997t iRos et alj|200qb. B07 1 2+47 2 
dJackson et all fl998l). B 1422+ 23 1 dPatnaik et alj [l999l [l992h 
B1933+503 dSvkes et al.lll998l) . and B2045+265 dMcKean et al.l 
l2007h PI 
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Fig. 1. The cumulative flux perturbation distribution for 5 ob- 
served lens systems. Solid lines represent positive parity images, 
while dashed represent negative parity images. The brighter of 
two images with the same parity is denoted by a thick line and 
the fainter by a thin line. 

Using our aut omated lens modeling cod e, we try to repro- 
duce the results of lKochanek & Da lai (2004). As discussed pre- 



2 We fit in the source plane by minimizing the dispersion in the mod- 
eled source positions. Further details are found in AppendixlAl 

3 Additional data taken from http://cfa-www.harvard.edu/castles 
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viously, we fit the positions of the images and the lensing galaxy 
(10 constraints) with a SIE lens plus external shear. The fluxes 
(and magnifications) of the images are not used as constraints. 
The estimated flux perturbation for each image is defined by 



Table 1. Best-Fit Parameters for Observational Sample 



8f,i = log(/ob S; ,y/mod,/), 



(3) 



where / b s ,, is the observed flux of the z'th image and / mo d ; , is the 
modeled flux of the 2th image. The modeled fluxes are related to 
the modeled magnifications by 



fmod,i — S () //mod,/ 



(4) 



where So is the unlensed source flux. The unlensed source flux, 
however, is unknown. If there were no flux perturbations, the 
source flux would be equivalent to So = /obs,;/A ( mod,;- It can be 
estimated by 



1 2 



fohsj 



2 *r~^ j-^modj 



(5) 



where - in order to minimize the biasing due to flux perturba- 
tions - we discard the images with the largest and smallest ratios 
of /obs,;/yUmod,; and average over the remaining two images Q 

We separate the four images of each lens by brightness 
and parity, denoting positive and negative parity images by "+" 
and "— " signs and brightest and faintest images by and 
"*" signs. Thus, each lens has a brightest positive parity im- 
age (!"+■*•), a faintest positive parity image (J+«), a brightest 
negative parity image (/-*), and a faintest negative parity im- 
age (/-*). Correspondingly, each image, J;, has a magnifica- 
tion, Hi, and flux perturbation, <S/,,, or magnification perturbation, 
5i = logOu b S) i/jU mo d,i)- 

The cumulative distribution of flux perturbations, N(< 
(!>/_;) /(total number of lenses), is shown in Figure Q] The results 
are plotted so that negative (positive) values represent images 
with fluxes that are dimmer (brighter) than the best-fit modeled 
fluxes. The results do no t reproduce the exact relation seen in 
Kochane k & Dalall d2004l) : specifically, our results show overall 
larger magnification perturbations in the faint images (in particu- 
lar in J-*). In ad dition, not all of our lenses show that 5/-* < 
as is the case in iKochanek & Dalall (|2004). This does not con- 
stitute a discrepancy since the unlensed source flux has been 
estimated slightly differe ntly. Our results generally agree with 
IKochanek & Dalall d2004l) as parity dependence is very clearly 
seen and is such that 5/,-* < £>/, ail others (i-e., the brightest neg- 
ative parity image is less magnified or more demagnified rel- 
ative to its modeled magnification than other images are com- 
pared to their modeled magnifications). Our c ode uses com- 
monly used methods for lens modeling (see, e.g., Keeton 2001a) 
and commonly used lens paramete r izatio ns but is independent of 
the analysis by Kochanek & Dalai ( 2004), suggesting that parity 
dependence is a generic result of parametric lens modeling. In 
addition, the errors in the fluxes (from less than one percent to 
10%) are not large enough to obliterate the observed effect. 

The best-fit parameters for b, q, y, and 6 y are shown in Table 
Q] We use these parameters, in addition to the halo orientation 
and the best-fit source and lens positions to create a set of smooth 
models on which to test substructure models. 



Lens 


b (") 


q 


y 


By (rad) " 


MG0414+0534 


1.37 


0.75 


0.11 


0.98 


B0712+472 


1.05 


0.46 


0.10 


1.73 


B 1422+231 


0.94 


0.69 


0.16 


0.07 


B 1933+503 


0.67 


0.60 


0.06 


-0.22 


B2045+265 


1.36 


0.63 


0.06 


0.62 



" Shear orientation is relative to halo main axis. 

4. Mock Lenses 

4.1. Generating Mock Observations 

We test substructure models with mock observations of lens sys- 
tems. We begin with a set of models for lens halos and lens envi- 
ronments: the best-fit smooth macromodels of the previous sec- 
tion as summarized in Table Q] To the lens halos we add clumps 
as discussed in the following section. 

For a single mock observation of a lens macromodel and sub- 
structure, we choose a source position and find the associated 
image positions and magnifications by solving the lens equation 
using a Newton-Raphson method on a grid of possible image 
positions. Gaussian observational errors are added to image po- 
sitions and lens galaxy position: we adopt an observational error 
of 3 mas. Systems that result in 4 "observed" images are mod- 
eled using our automated lensing code, and the results are used 
to create a cumulative magnification perturbation distribution. 

When we choose a source position, we have two options. We 
can use the best-fit source positions of the observational sample. 
This results in a lens configuration (e.g., cusp, fold, or cross) that 
is similar to that of the observed lenses. Alternately, we can sam- 
ple the source plane and create a variety of lens configurations. 
When sampling the source plane, we account for magnification 
bias, assigning source po sitions by sampling t he im age plane 
uniformly as described bv lKeeton & Zabludoffl (i2004)rl Results 
where the source plane is sampled with a uniform weighting in 
the image plane are labeled 'UW.' Regardless of how the source 
positions are assigned, we create and model -1000 realizations 
for each lens macromodel, for a total of ~5000 realizations. 

4.2. Substructure Models 

In order to approximate the clumpy distributions of matter we 
expect to find in galaxy halos, we add dark matter substruc- 
ture to the smooth macromodel. These clumps are modeled as 
projected Moore-like profiles (see, for comparison, Moore et al. 
1999) with density profile 



k(x) 



K S (X 112 + X 2 )- 





1 ,R<n 

,R > r, 



(6) 



where x — R/r s , R is the projected separation from the center of 
the clump, r s is the scale radius, and r t is the tidal radius of the 
clump. We test substructure models based on simulations and ad 
hoc substructure models as described in the following sections. 

4.2.1 . iMao et al.l j2004|) Substructure Model 

Lensing is most s ensitive to the clu mps near the Einstein radius 
of the lens galaxv. lMao et al. I (l2004l) estimate that the Einstein ra- 
dius is ~3% of the virial radius of a lensing halo (b = 0.03/? v ir)- 



IKochanek & D alai ( 2004) average over all four images. 



5 iKeeton & Zabludoff (2004) show that for sources with a power law 
luminosity function, dN/dS oc 5 ~ v with v = 2, magnification weighting 
in the source plane is equivalent to uniform sampling of the image plane. 
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If we assume that the halos used in the simulations of Mao et ail 
(2004) can be modeled by SIEs, then the mass within b is also 
~3% of the halo mass (M(< b) - 0.03M ha i o )- Given a surface 
mass density in subclumps at the Einstein radius of ~0.5%, the 
total mass in subclu mps is ~ 10~ 4 M hai o . The mass resolution of 
the simulations of iMao et al.l d2004l) is also 10" 4 Mh a i o , and sim- 
ulations have shown that smaller clumps are significantly more 
comm on than larger clu mps. We can then assume that the re- 
sult o f lMao et all (2004) can be approximated by a substructure 
model where 1 clump of 10~ 4 Mh a i o is centered on the Einstein 
radius. 

IMao et ail ([2004) also find that most clumps have a tidal ra- 
dius which is ~0. 1 of the halo virial radius. Thus, clumps that are 
found at separations up to four times the Einstein radius, 4b, con- 
tribute to the measurable substructure mass density. A better sub- 
structure model would inc lude subclumps w ithin 4b of the lens 
center. Within that radius, IMao et ai] (|2004) find a substructure 
mass fraction of ~2%. Our fiducial substructure model, then, in- 
cludes 24 subclumps of 10~ 4 Mh a i o within 4b of the lens center. 
Each clump has a tidal radius equal to 0. 1 of the halo virial ra- 
dius. Fiducially, we set the scale radius to b/20. The clumps are 
placed to follow the SIE density distribution. We refer to this 
model subsequently as M04. 



The cumulative distribution of magnification perturbations, 
N(< Si)/N tot , using the fiducial M04 model are plotted in Figure 
2] Here, w e can s ee results that appear to be consistent with 
Rozo et al.l d2006l) : the median perturbation is near zero, al- 
though the tails are not negligible. 
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4.2.2. The 'ALL Substructure Model 

As the substructure model based on numerical simulations will 
prove insufficient to reproduce the observed parity dependence, 
we create a substructure model with greater effect on the magni- 
fication perturbation. This model places one clump within a scale 
radius of each of the positions where images would be found 
in the absence of substructure. We refer to this substructure 
model as ALL' since all images have a corresponding clump. 
Fiducially, we use the same tidal radius, mass, and scale radius 
as in the M04 model. 

5. Results 

When we model real observations in Section [3j we are obliged 
to compare observed fluxes to modeled fluxes. In our mock ob- 
servations, however, we can 'observe' magnifications and, thus, 
we compare observed and modeled magnifications. The magni- 
fication perturbation is estimated by 

5i = \og(fi ohs j/fi modJ ), (7) 

where yU b s ,, is the observed magnification of the ith image and 
fJ-modj is the modeled magnification of the ith image 

We estimate the magnification perturbations for three differ- 
ent scenarios for substructure in the following section. First, we 
see if substructure models based upon dark matter simulations 
are sufficient to induce parity dependence. Second, we investi- 
gate whether satellite galaxies in lens halos can account for the 
observed parity dependence. Finally, we design a substructure 
model that does induce parity dependence and discuss the impli- 
cations of such a model. 

5.1. Comparison to \Mao etal\ \2004) Substructure Model 

We investi gate if the s ubstru cture constraints found in the sim- 
ulations of Ma o et all (|2004|) can reproduce the observed parity 
dependence in this section. We first test the fiducial M04 model 
described in Section 14.2.11 Then, explore the effects of vary- 
ing substructure model parameters individually: r s , r t , subclump 
mass, and / su b. 
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Fig. 2. Parity dependence in the cumulative distribution of mag- 
nification pertu rbations using the substructure model following 
IMao et ail(l2004l) . The top panel shows the observed lens config- 
urations while the bottom panel varies the source positions used 
(see text). Solid lines represent positive parity images, while 
dashed represent negative parity images. The brighter of two 
images with the same parity is denoted by a thick line and the 
fainter by a thin line. 

Given the small size of the observational sample, the excep- 
tionality of the observed parity dependence is best quantified by 
estimating the probability of finding a sample of 5 lenses with 
significant parity dependence. Increasing the number of realiza- 
tions modeled such that each lens macromodel is used at least 
5000 times, we draw 5000 sets of 5 lenses (one mock observa- 
tion for each of the lens macromodels). We create test criteria for 
observing parity dependence in the sets. One important aspect in 
developing criteria is specifying how different the magnification 
perturbations of and I are in order to exclude sets where 
the difference is negligible and no parity dependence would be 
observed. (5/ i+ * - £>/,-★ = 0.06 is the smallest difference seen 
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in the observational sample, and we adopt this in the criteria for 
S + + - The three sets of criteria we test are: 

1 . <S_* < 6++ everywhere and at least 60% of the distribution 
shows 6 + + - i5_* > 0.06 

2. < <5 a n others everywhere and at least 60% of the distribu- 
tion shows 6 + + - <5_+ > 0.06 



Table 2. Results Using Substructure Models Based on Mao et al. 
(2004) 



6-+ < <5aii others and 6++ - 
dependence probability") 



6-+ > 0.06 everywhere ("parity 



Of the 5000 sets, 170 (3%) satisfy the first set of criteria, 58 sets 
(1%) satisfy the second, and 14 sets (0.3%) satisfy the third. So, 
given the M04 substructure model, the probability of finding a 
parity dependence like the one observed is small but maybe not 
insignificant. We designate the third set of criteria as the par- 
ity dependence probability. This measure is more useful than a 
standard Kolmogorov-Smirnov test, as we are most interested in 
the relation between magnification perturbations of different im- 
ages, not comparing the absolute distribution of observed flux 
perturbations to magnification perturbations of mock observa- 
tions. 

In the bottom panel of Figure [2] we vary the source posi- 
tions as described in Section 14.11 Here, the parity dependence 
probability is 5 out of 5000 sets (0.1%). The particular lens con- 
figurations in the observed sample, then, seem to be mildly more 
effective in creating parity dependence than a sampling of source 
positions. Henceforth, we use a smaller number of sets to calcu- 
late the probability, a total of 500. 
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Fig. 3. Varying the scale radius and tidal radius. The left-hand 
panels test larger values of the scale radius while the right-hand 
panels test smaller values of the tidal radius. The parity depen- 
dence probability for each test is in parentheses. Solid lines rep- 
resent positive parity images, while dashed represent negative 
parity images. The brighter of two images with the same parity 
is denoted by a thick line and the fainter by a thin line. 



We have cho sen to parameterize the substructure model of 
iMao et all d2004l) in a reasonable, but not unique, way. Here we 
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test how the parity dependence depends on the parameters cho- 
sen, and results are summarized in Table [2] We vary the tidal 
radius r, or the scale radius r s individually, while keeping the 
other subclump parameters fixed at fiducial values in Figure [3] 
Increasing the scale radius reduces the parity dependence in- 
duced by substructure, as the mass within the scale radius is 
spread out over a larger area. Decreasing the tidal radius puts 
more mass into a smaller area and increases the parity depen- 
dence induced by substructure. The parity dependence proba- 
bility increases by an order of magnitude to ~ 2%, but remains 
small. Varying the scale radius or the tidal radius does not change 
the parity dependence probability in a linear way. For example, 
the smallest tidal radius, lO -3 /?^, is smaller than the fiducial 
scale radius and has the same parity dependence probability as a 
substructure model with tidal radius of 10 2 /? v ir- 

We do an additional test, modeling the subclumps as SISs 
which has steeper inner profiles than our Moore-like profiles and 
which might better approximate clumps with satellite galaxies in 
their centers. At fixed mass, however, the perturbers are no better 
than the M04 in inducing parity dependence. 
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Fig. 4. Test using a realistic mass function. The parity depen- 
dence probability is in parentheses. Solid lines represent positive 
parity images, while dashed represent negative parity images. 
The brighter of two images with the same parity is denoted by a 
thick line and the fainter by a thin line. 
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We try a slightly more realistic version of the substructure 
model, including a mass function for su bclumps. Dark mat - 
ter simulations suggest that n(m) ~ mr x i dGhigna et alJ uOOO). 
The scale radius is scaled by the square root of the mass: at 
M = l(T 4 Mhaio, r s has the value used in the fiducial model, and 
it is larger for larger masses and smaller for smaller masses. The 
results are shown in Figure [4] where it is clear that this more re- 
alistic model is no more effective in inducing parity dependence 
than the fiducial model. 
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Fig. 5. Testing a larger substructure fraction (/ su b = 25% within 
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clumps (bottom). The parity dependence probability for each test 
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substructure distribution at z = 0. Since a significant number of 
strong lenses are found nearer to z = 1, a significantly larger 
amount of substructure is expected. If we increase the substruc- 
ture mass fraction, the size of the magnification perturbations 
correspondingly increases, as does the probability of observed 
parity dependence. This can be seen in Figure [5] where 300 
clumps are used (which are 12.5 times more than in the fidu- 
cial case and correspond to a substructure mass fraction within 
4b of 25%). Using a larger substructure mass fraction, we also 
test the effect of substructure mass. Keeping both the tidal radius 
and scale radius fixed, fewer, larger mass clumps are more effec- 
tive than a greater number of smaller mass clumps, but the parity 
dependence probability remains no greater than a few percent. 

Substructure models based on dark matter simulations, then, 
seem insufficient to account for the observed parity dependence 
unless the observed parity dependence is an outlier in the distri- 
bution. If, however, the substructure mass fraction is an order of 
magnitude larger, the probability of the observed parity depen- 
dence may be a few percent. In addition, clumps that are well- 
concentrated have a greater effect than those that are less concen- 
trated. Well-concentrated clump might be found in higher reso- 
lution dark matter simulations: clumps that are below the mass 
resolution and/or are tidally destroyed faster in lower resolution 
simulations. In addition, in simulations that include the effects of 
baryons, we may see longer surviving subclumps that have their 
central densities boosted by cold baryonic components. 



5.2. The Effect of Luminous Satellites 

We find it difficult to reproduce the observed parity dependence 
using realistic models of dark matter substructure. However, 
of the 5 lenses in our sample, 2 of them (MG04 14+05 34, 
B2045+503) have luminous companions that lie in and around 
the Einstein radius of the lens halo. These companions may be 
satellite galaxies in the lens plane0 We test the effect of placing 
a single luminous satellite within two Einstein radii of the lens 
center in Figure[6] The satellite is modeled as a SIS with trunca- 
tion radius set to the Einstein radius of the lens halo (r tiSa t = b). 
The Einstein radius of the SIS is one-tenth the Einstein radius 
of the lens (b sat = b/10) and it has a mass of 0.01Mh a i o - This is 
similar to the best-fit lens models for the companion objects in 
MG0 414+0534 and B2045+503 dRos et al.ll2000tlMcKean et al.1 
12007b which parameterize the companions as untruncated SISs 
and have best-fit values of b sat = 0.17 and 0.08, respectively. In 
our case, while the luminous satellite seems to increase the ab- 
solute size of the magnification perturbations, the probability of 
observing parity dependence is estimated to be as low as in the 
M04 substructure model (1 out of 500 sets of realizations). 

If the additional mass were much larger, it could 1 .) steepen 
the inner profile of the lens from a SIE and 2.) significantly in- 
crease the Einstein radius of the lens. This would induce a parity 
dependence such that all <5, < and < S + + . However, faint 
negative parity images would be just as or more demagnified 
compared to bright negative parity images. 

In a statistical sense, luminous satellites do not appear to 
resolve the observed parity dependence. Regardless, we addi- 
tionally check to see if they could eliminate the observed par- 



The amount of substructure in the M04 models seem insuf- 
ficient to match the observational results, unless our observa- 
tional results represent an outlier in the distribution of magni- 
fication perturbations. The M04 model, however, represents the 



6 The prevalence of luminous satellites in anomalous flux ratio 
systems has been noted previously: 5 of the 22 gravitat ional lenses 
found in the Cosmic Le ns All-Sky Survey (CLASS; iBrowne et alJ 
120031 : iMvers et all l2003h have luminous companions (B 1608+656, 
B2045+503, MG0414+0534, B1127+385, B1359+154). 
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Fig. 6. Testing the effect of a single luminous satellite. The par- 
ity dependence probability is in parentheses. Solid lines repre- 
sent positive parity images, while dashed represent negative par- 
ity images. The brighter of two images with the same parity is 
denoted by a thick line and the fainter by a thin line. 

ity dependence in the specific cases of MG0414+0534 and 
B2045+503. 



all best-fit modeled magnifications find that fi_+ > . The 
brightest images in the system are the fold images, and, given 
the distance of the satellite object, it is not surprising that it does 
not seem to effect the ratio of those fluxes. 

B2045+503 exhibits a cusp configuration with three im- 
ages positioned near each other. Its source and lens redshifts 
have been measured at z s = 1.28 and z\ = 0.867, re- 
spectively dFassnacht et al.| [l999). The companion galaxy (G2) 
lies between the lensing galaxy (Gl) and the cusp images 
and has a photometric reds hift consistent with the lensing 
galaxy (|McKean et al.l 120071) . The best-fit smooth model by 
iMcKean et all (120071) for Gl and G2 finds that the ellipticity 
of Gl is oriented 90 degrees away from the cusp images while 
a large external shear is pointed along the line from Gl to the 
cusp images. In addition, their best-fit model for G2 is very el- 
liptical (q - 0.133) and pointed along the same direction as the 
shear. Including the previously found best-fit parameters for G2 
- Einstein radius, axis ratio, orientation, and position - we can 
significantly reduce the x 1 of our original best-fit, but the mag- 
nification perturbations are not reduced in this case: <5/,-* < 
and < £/, a ii others- As in the case of MG0414+034, the 

observed fluxes show that I + + is brighter than J_„, while the 
best-fit modeled magnifi cations find that ^u_+ > We use 

the ob served fluxes from Mc Kean et al.l (l2007l) .lKoopmans et al.l 
(120031) find that B2045+503 is affected by variability in the radio 
fluxes, however, our results are relatively unaffected. 

Luminous satellites, then, do not appear sufficient to explain 
the observed parity dependence either in the statistical sense or 
in the particular cases in our observational sample. 



Table 3. Flux Perturbations in MG0414+0534 and B2045+503 



Lens and Model 


s f,+* 


<5/, + . 




6f,-, 


MG04 14+0534 


-0.09 


0.11 


-0.16 


0.07 


include sat 


0.013 


0.12 


-0.052 


-0.013 


B2045+503 


-0.18 


-0.12 


-0.53 


0.57 


include sat 


-0.27 


0.16 


-0.73 


0.70 



MG0414+0534 exhibits a fold configuration - two of the 
four images are found close to each other. The source and lens 
redshifts are found at z s = 2.6 and zi = 0-9 6, respectively 
dLawrence et al. I ll 9951: iTonrv & Kochanek| [r999). The compan- 
ion object lies near the Einstein radius and far from the fold im- 
ages. It does not have a measured redshift, but we assume that 
it lies in the lens plane. We find an initial smooth lens model, 
using 10 constraints and 9 parameters. We add the companion 
object to the lens model as a SIS, fixing its position to match 
observations and its Einstein radius to 0.12 times the Einstein 
radius of the lensQ We add no new parameters nor constraints 
and refit, finding a new, smaller^- 2 . The magnification perturba- 
tions for each image are shown in Table [3] Here the inclusion 
of the satellite object clearly reduces the size of the magnifi- 
cation perturbations. While in the original fit, all magnification 
perturbations were ~25%, now two of the four are less than 5%. 
However, the size of 6 1++ - 6 / remains consistent. In addi- 
tion, 6f-.+ < (5/, a ii others- This parity dependence appears robust: 
all observed fluxes show that 1 is brighter than J_„ while 

7 Ros et al. (2000) find a larger value of the Einstein radius of the 
satellite than we use. We find that larger values for the Einstein radius 
induce extra ima ges in the lens. Note also that we us e the fluxes from 
iKatz et al.1 dl997i) and not those in the dRos et alJi OOO). although results 
are consistent. 



5.3. Creating an Effective Substructure Model ('ALL') 

Given the difficulty in reproducing the observed parity depen- 
dence with well-motivated models of dark and luminous sub- 
structure, a pertinent question to ask is what substructure model 
will induce parity dependence. There exists a subsample of M04 
substructure realizations which do show appreciable parity de- 
pendence - the fraction of those realizations with at least one 
clump located within a scale radius of an observed image. This 
sample has a parity dependence probability of ~14%; however, 
the fraction of realizations in which the M04 substructure model 
results in a close alignment of clump and image is small, ~ 4%. 

Based on that result, we create a substructure model with 
greater effect on the magnification perturbation. This model 
places one clump within a scale radius of each of the image po- 
sitions that would be found in the absence of substructure. As 
previously noted, we refer to this substructure model as 'ALL 
because all images have a corresponding clump. 

The results featuring a clear parity dependence are shown in 
Fi gure[7] Qualitatively , the di stributions are similar to the results 
of iKochanek & Dalall d2004l) : 6-+ < <5 + * while the distribu- 
tion of 5 + + and 6f a i nt images are very similar. The UW test has 
similar results although the faint images do not correspond to 
6++ as in the previous case. The probability of observing parity 
dependence has increased significantly to 40% and 80% when 
sampling the source plane and when using specific source po- 
sitions, respectively. A similar model - placing small, truncated 
SISs near images - has been tested by Kochanek & Dalai (2004) 
with similar qualitative results. 

The efficacy of this model suggest that the requirements 
for parity dependence may be significant numbers of well- 
concentrated and well-aligned perturbers. It also supports the 
possibility that m ore substructures than the amount found by 
iMao et all d2004l) are required to achieve the observed parity de- 
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Fig. 7. Parity dependence when each image in a lens is paired 
with a nearby subclump. The top panel shows results using the 
best-fit source positions from the original data, while the bottom 
panel samples the source plane. The parity dependence probabil- 
ity is in parentheses. Solid lines represent positive parity images, 
while dashed represent negative parity images. The brighter of 
two images with the same parity is denoted by a thick line and 
the fainter by a thin line. 
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Fig. 8. Testing the effect of individual clumps on the parity de- 
pendence. For the top 4 panels, clockwise from the top-left, 
each panel shows the results for placing one clump only near 
the bright positive parity image (BP), faint positive parity image 
(FP), faint negative parity images (FN) and bright negative par- 
ity image (BN). The bottom-right panel shows the effect of plac- 
ing one clump near a random image position and the bottom-left 
panel shows the effect of placing two clumps near two random 
image positions. The parity dependence probability is in paren- 
theses. Solid lines represent positive parity images, while dashed 
represent negative parity images. The brighter of two images 
with the same parity is denoted by a thick line and the fainter 
by a thin line. 



pendence. While the ALL substructure model uses fewer sub- 
clumps and, therefore, ha s a smaller meas urable substructure 
mass fraction (l/6th of the iMao et al.l (2004) result), the place- 
ment of clumps is perverse. We expect that, in a realistic sub- 
structure model, the placement of clumps in the projected plane 
is stochastic, and the probability of a chance alignment between 
image and clump is small. 

Now that we have a substructure model that induces the ob- 
served parity dependence, we test which of the clumps are the 
most effective in Figure [8] Clumps near faint images have rel- 
atively mild effects; clumps near bright images show larger ef- 
fects on the magnification perturbation distribution. In particular, 
clumps near 1 cause a significant demagnification in that im- 
age, while clumps near I + + cause large negative magnification 
perturbations in that image, as well as a pronounced large de- 
magnification tail in J_*. Intriguingly, although at least some 
of the lens configurations have and I as a close pair, 
clumps near I_+ have little effect on while the converse 



is not true for clumps near In the bottom panels of Fig. [8] 
one chance alignment between image and clump is sufficient to 
boost the probability of observing parity dependence, and two 
chance alignments puts the probability at nearly 50%. 

The results of this model can be described simply by inspect- 
ing the eigenvalues of the magnification matrix. The lens equa- 
tion describes the relation between the source and image posi- 
tions u = x - Viff(x), where u is the source position, x is the 
image position, and if/ is the lensing potential. The magnification 
matrix is given by du/dx. The eigenvalues of the magnification 
matrix are given by 

M = l-K-y (8) 

and 

X 2 = 1 - k + y, (9) 
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and the magnification is // = l/(AiA2). In the case of double- 
positive parity images, both A\ > and Ai > 0, while in the 
case of negative parity images Ai < and A2 > 0. In the case 
of the brightest images - for example, a fold - we would expect 
that A\ ~ for both images. By placing a clump near the im- 
age, we add a little to the surface mass density k. For the positive 
parity image, in general, this makes A\ smaller and closer to zero 
and increases the magnification. For the negative parity image, in 
general, this makes A\ more negative and further from zero, de- 
creasing the magnification. Thus, generically, adding mass near 
images should make 1.) the brightest negative parity image dim- 
mer than a smooth model would suggest (£_.*. < 0) and 2.) make 
the brightest positive parity image brighter than a smooth model 
would suggest > 0). 




40 60 



100 



Fig. 9. The cumulative distribution of the total magnification for 
realizations using the ALL substructure model and sampling the 
source positions (dashed line). The solid line shows the distribu- 
tion of total magnification using no substructure. 



As previously noted, a substructure model such as ALL is not 
plausible. But it might be realistic if it induced a large increase 
in the lensing probability. Two possibilities are that clumps have 
increased the size of the tangential caustic. Given that the mass 
fraction in substructures added by the ALL model is very small, 
this possibility is unlikely. The other possibility is that the total 
magnification of the lensed systems is increased making the lens 
system more likely to be observed. This seems highly unlikely 
given that the main effect is to demagnify one image greatly. 
Nonetheless we test this possibility using the ALL substructure 
model and sampling the source positions in Figure [9] where we 
plot the cumulative distribution of total magnification. At the 
smallest magnifications, the ALL substructure model has larger 
total magnification than in the no substructure case. This is prob- 
ably the result of the lenses having more mass in the ALL case 
than in the no substructure case and does not reflect any real 
magnification bias. The cumulative distribution in the ALL case 
rises more steeply and, overall, the result is that substructure ac- 
tually decreases the total magnification. 



6. Discussion 

We discuss here the different contributions to the observed parity 
dependence from 1.) lens modeling and 2.) astrometric perturba- 
tions. 

In general, we find lens modeling finds appropriate best-fit 
parameters (see Figure [Tol l. The distribution of best-fitparam- 
eters peaks around the no-substructure input parameter^ In ad- 
dition, the width of the best-fit parameter values from the M04 
and from the ALL substructure models are similar. The widths 
of the distribution of best-fits, however, can be fairly large (10- 
20%). A bias in the lens modeling code does not appear to be a 
contributing factor to the parity dependence. 
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Fig. 10. The best-fit smooth model parameters using M04 (solid 
lines) and ALL substructure models (dashed lines) in conjunc- 
tion with the best-fit macromodel for MG0414+0534 (see Table 
[T). Source points are sampled. The input macromodel values are 
shown as vertical lines. 



In fact, lens modeling as a whole has only a small effect on 
the parity dependence. In FigureQT] we perform two tests on the 
magnifications using no lens modeling and sampling the source 
positions. In the left-hand panel, we show the distribution of 
magnification perturbations using the ALL substructure model 
with no lens modeling at all. Instead of comparing the observed 
and modeled image magnifications, we compare the magnifica- 
tion at the image position of a mock lens with substructure with 
the magnification at the image position of a similar mock lens 
but without substructure. This result includes the effect of shifts 
in the image positions due to substructure: 



(W shift = log 



/^obs,sub(-^obs,sub) 
/^obs,no sub(-^obs,no sub) 



(10) 



8 The distributions for the Einstein radius b peak slightly larger than 
the input value, which is not unexpected given that we added mass to 
the lenses in the form of clumps. 
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Fig. 11. Comparing the magnification of images with and with- 
out substructure (no lens modeling). In the substructure case, 
each image in a lens is paired with a nearby subclump. Left: The 
magnification at the image position in a clumpy lens compared to 
the magnification at the image position of a smooth lens. Right: 
The magnification at the image position in a clumpy lens com- 
pared to the magnification at the same image position but using a 
smooth lens model. Solid lines represent positive parity images, 
while dashed represent negative parity images. The brighter of 
two images with the same parity is denoted by a thick line and 
the fainter by a thin line. 

In the right-hand panel, we compare the magnification at the im- 
age position of a mock lens with substructure with the calculated 
magnification at the same image position (the image position 
given substructure) but putting no substructure into the mock 
lens. This result, then, has no shift in image position. This test 
is similar to simple explanation for parity dependence that we 
described in Section |531 



shift = log 



/^obs,sub(-^obs,sub) 
/^obs.no sub(-^obs,sub) 



(ID 



In both cases, strong parity dependence of the kind seen in Fig. 
|7]is evident even though we have done no lensing modeling. 
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Fig. 12. Measuring the astrometric perturbations with no lens 
modeling.. Left: The results for placing one clump near each im- 
age (ALL)^ip/j£^The results using a substructure model based 
on iMao et al.l (120041) (M04). Solid lines represent positive par- 
ity images, while dashed represent negative parity images. The 
brighter of two images with the same parity is denoted by a thick 
line and the fainter by a thin line. 

We examine the size of the shifts in image position due to 
substructure, and without any lens modeling, in Figure [12] The 
differences between using the M04 substructure model and the 
ALL substructure model seem modest. In both cases, the median 



change in image position between a lens system with no sub- 
structure and one with substructure is ~30 mas. However, when 
examining the images separately, we see a significant difference 
between substructure models. For the ALL substructure model, 
the astrometric perturbations for are significantly smaller 
than that of In the M04 case, these two cases are similarly 
large. 

As mentioned previously, we add clumps to the positions 
where the input macromodel would find images in the absence of 
substructure. From Figure [12] we see that whether the observed 
image position is near a clump depends on the image. For ex- 
ample, in the ALL model, it seems that a significant number of 
faint images and of may remain well aligned with their per- 
turbers, while the observed I + + may have moved a significant 
distance relative to the scale radius of the clump. 

7. Conclusions 

Radio observations of quasar sources strongly lensed by galaxy 
lenses provide information on the mass distribution along the 
line-of-sight and, as such, have been suggested as probes of 
galaxy substructure and a small-scale test of CDM. One sim- 
ple measure of the subclumps is the substructure mass fraction 
at the scale of the Einstein radius in a galaxy lens. Simulations 
have predicted that this fraction is / su b ~ 0.5% at z — 0. 

Previous comparisons of observations to simulations have 
suggested that there may be more substructure than pre- 
dicted in simulati o ns w it h varying le v els o f discrepancy 
dDalal & Kochanekl l2002t iMaccio et all 120061: lAmara et al.l 
120061) . I Kochanek & Da lai (2004) pointed out that such magnifi- 
cation perturbations did not behave identically in each image and 
the behavior depended upon the parity (orientation) of the im- 
age. The observed parity dependence they found is such that the 
brightest negative parity image is demagnified with respect to the 
magnification predicted by smooth models and has a magnifica- 
tion perturbation smaller than that of the other images (5_* < 
and <5_+ < £ a ii others)- Here, we use the observed parity depen- 
dence in order to extract additional information about the sub- 
structure distribution. 

We create mock observations of strong lens systems with 
clumpy lens halos and compare the distribution of magnification 
perturbations between the mock and observed samples using an 
automated lens modeling code. Overall, our results are consis- 
tent with previous results which suggest that there is more sub- 
structure in the observations than can be accounted for in dark 
matter simulations. Our conclusions are summarized as follows: 

1. Parity dependence in a sample of four-image lenses 
is found such that the brightest negative parity has a 
magnification perturbation smaller (less magnified or more 
demagnified) than that of the other images ( cL* < 5aii others)- 

2. This effect is difficult to reproduce using substructure 
models based on dark matter simulations; the probability 
of generating a sample with properties like the observed 
sample is less than 1%. 

3. While a significant fraction of lens galaxies have been found 
to have luminous companions, satellite galaxies do not 
appear sufficient to induce parity dependence. 

4. Parity dependence requires that images and clumps are 
well-aligned and clumps are sufficiently concentrated. 
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5. More substructure than predicted in simulations may be re- 
quired, either by resolving smaller but well-concentrated 
clumps or by increasing the overall substructure mass frac- 
tion. 

Given the small sample of observed radio lenses, it is diffi- 
cult to put robust constraints on the size or significance of the 
observed parity dependence. In addition, while we approximate 
the observed parity dependence using simple substructure mod- 
els, we cannot exclude more complicated models for the lens 
halos. We also consider only substructure within the lens halo 
itself and do not consider the possible parity dependence effects 
of the significant numbers of clumps that may be found along 
the line-of-sight. Previous studies of the imp ortance of line-of- 
sight clumps h ave found conflicting results dChen et al.| [2003; 
Metcalf 2005a b). The large amount of line-of-sight structure 
suggests that this may be an effective way of resolving the parity 
dependence problem. However, line-of-sight clumps are less ef- 
fective than substructure in inducing magnification perturbations 
and the overall effect of line-of-sight clumps is unclear. 

Parity dependence as observed in four-image lenses with ra- 
dio observations is a strong measure of anomalous fluxes and a 
constraint on substructure in galaxies. In this paper, we describe 
the contradictions between the observed parity dependence and 
previous theoretical estimates and resolve the discrepancies us- 
ing modeling of mock observations of lenses. We attempt to 
constrain the kinds of substructure models that will reproduce 
the observed parity dependence. The substructure models nec- 
essary to approximate the observed results is consistent with 
suggestions that there is more substructure in galaxies than in 
estimates from simulations. In addition, we show that it is un- 
likely that anomalous flux ratios are due solely to observed lu- 
minous satellite galaxies. Future studies of CDM substructure 
using strong gravitational lenses may require larger samples of 
observed lenses as well as more physically realistic simulations 
of dark matter and baryons. 
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Appendix A: Lens Modeling Algorithm 

Our lens modeling algorithm consists of 7 steps (some of which 
are repeated): 

1. Fix the value for the Einstein radius to half the maximum 
separation of images. Set y, 6 y to zero. Test values of q and 
Q q on a grid in the source plane. 

2. Using the best fit value of q and 9 q , fit for b, q, and 6 q in the 
source plane. 

3. Using the best fit values of b, q and 9 q , test values of y and 
6 y on a grid in the source plane. 

4. Fit for b, q, 6 q , y, and 6 y in the source plane. 

5. Do three iterations of fitting b, q, 9 q , y, y , x cente r, and 
Reenter in the source plane. Each trial perturbs the image 
positions slightly. The set of parameters with the smallest 
X 1 are accepted as the initial guess parameters for the next 
iteration. If the x 2 increases, those best-fit parameters are 



rejected. 

6. Do up to 10 iterations of fitting in the image plane. In each 
iteration, the image positions are perturbed slightly. The set 
of parameters with the smallest x 2 are accepted as the initial 
guess parameters for the next iteration. If the x 2 increases, 
those best-fit parameters are rejected. Iterations are stopped 
if x 2 IN, the chi-squared per number of constraints, drops 
below 1. 

7. Check the x 1 f° r the best-fit parameters using original, ob- 
served image positions. Recheck modeled image positions 
for best-fit parameters. If the wrong number of images or 
parity violation is found, start again (up to 6 times) and set 
different values for b and q in Step 3. 



The source plane x 2 is calculated as in iKeetonl (12001 a) and 
reproduced here for convenience: 



where 



R 



o-l 



Ri 



6U = Hnbs - H„ 



COS 0, T 



(A.l) 



(A.2) 



(A.3) 



(A.4) 



H is the magnification matrix and u b s ,i is the source position 
counterpoint to each observed image. 

Fitting in the source plane avoids the need to solve the lens 
equation and allows the best model source position to be found 
analytically: 



"mod — A ■ b, 

where 

i 

and 



7 ST 1 



The image plane fitting uses &x 



x 2 

A img 



= Yl 6X J • S i 1 ■ 6Xi + Center ' S'i 



■ 5x c 



where 
and 

(^center,/ — X cen t ero bs — X cen t e r im od 



(A.5) 
(A.6) 

(A.7) 

(A.8) 

(A.9) 
(A. 10) 



and the sum extends over all images, x bs,; and x mo d,,- are the 
observed and modeled positions of image i, and x center>0 b s and 
Xcenter.mod are the observed and modeled positions of the center 
of the lens potential. The astrometric uncertainties for image i 
are described by the covariance matrix. The covariance matrix 
we use is 



erf 
aj 



(A.ll) 
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